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I. INTRODUCTION 

In 1982,JUtshuler, Aronov, and Khmelnitsky (AAK) 
established^ that electron-electron scattering in met- 
als is characterized by three (generally, distinct) time 
scales. These scales are phase-relaxation time t^, energy- 
relaxation time te, and out-scattering time r e . The for- 
mer one is quantum-mechanical and has no classical ana- 
log, while the two latter have a semi-classical interpreta- 
tion in terms of Boltzmann equation. The three scales 
differ in the case when the energy transferred between 
electrons in one collision is small as compared to the tem- 
perature of the system T. 

One can understand the difference between te and r e 
by considering the inelastic collision integral in Boltz- 
mann equation, 



St {/(e)} = J dexdujK{u)) 



(1) 



x { - /(e) [1 - /(e - uj)] /( £l ) [1 - /( Cl + w)] (out) 

+ [1 - /(e)] /(e - uj) [1 - fie,)] f(e 1 + w)}, (in) 

where /(e) is the electron distribution function, and the 
kernel K(uj) characterizes matrix elements of the inter- 
action, with the energy transfer uj. In clean 2D and 3-D 
systems, K(u>) is independent on the transmitted energy 
uj, K(tu) ~ 1/eF- This results in a Fermi liquid behavior 
of the inelastic rate l/r^ ~ max(e, T) 2 /ej?. The .situa- 



tion in disordered systems, however, is different^! u, the 
kernel K(u>) grows with the decrease of the transmitted 
frequency uj, 



K(uj) 
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where g(L u ) is the dimensionless conductance (in units 
if e 2 /-xTi) of the d - dimensional disordered sample of the 
size Lu = (D/uj) 1 / 2 , where D is the diffusion constant of 
the metallic sample. For more details on origin of Eq. (||), 
see e.g. Refs. ||7|. n 
Substituting Eq. (|) into Eq. (§), one cstimatesa 
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where e* = max(e, T). Equations (g) and (|3|) are appli- 
cable for systems in the metallic regime, g(L) 3> 1. In 
this regime e*TE 3> h~, i.e. quasiparticles are well-defined. 
Notice that, even though the kernel K is divergent, the 
energy relaxation rate (0) is finite because of the two en- 
ergy integrations in Eq. (|l|). Therefore, for the study of 
the phenomena governed by the Boltzmann equation, the 
infrared divergence of the matrix elements (|2|) does not 
cause any problems. These phenomena include, for in- 
stance, electron distribution function measured via tun- 
neling spectroscopy! ox crossover from 1/3 to y/3/4 shot 
noise in metallic wiresQ. 

It is not the end of the story, though. If we esti- 
mate only one ("out") term from the collision integral 
( l|) , we encounter an infrared divergence in two- and one- 
dimensional cases, 



SW{/(e)}' 



Ti 



T 



g{L T ) 



duj (T 

UJ \ UJ 



(4) 



where uj* is the low energy cut-off to be found, and 
Lt = \jDjT is the temperature length. (The same re- 
sult may be obtained from the calculation of the first loop 
correction to the self-energyQ.) This divergence of only 
one contribution to the collision integral is a simple con- 
sequence of the fact that each term in collision integral 
is not a gauge invariant quantity, and only both terms 
taken together have a physical meaning (|j|) , which is not 
cut-off dependent. One can argue, however, that r e has 
its own observable consequences for the quantum interfer- 
ence processes. Indeed, naive argument is that the "out" 
processes completely suppress the interference, whereas 
"in" processes are incalierent. Inclusion of some of the 
higher order processesQtl cures the divergence and makes 
the expression for l/r e finite. One may naively expect 
that T e found from such procedure is, indeed, responsible 
for the temperature behavior of quantum corrections. 

AAK showedE] that it is not correct for the tempera- 
ture behavior of weak localization correction, because the 
inelastic excitations with energy transfer smaller than de- 
coherence rate itself do not suppress this correction, see 
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Section 2.2.2 of Ref. and our Section |T| for the corre- 
sponding physical argument. This leads to the infrared 
cut-off ui* ~ l/r^ in Eq. (0) and to the self-consistency 
equation for the dephasing rate, 



T 4> 



g( L <t>Y 
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However, there is a prejudice, see e.g. Ref. 0, that the 
inelastic time governing the magnitude of the conduc- 
tance fluctuation is given by r e -C r^, so that r e has its 
own observable effect. 

In this paper we revisit this problem. We will show 
that the inelastic rate governing the mesoscopic fluctua- 
tions is precisely the same as for the weak localization, 
see Eq. (0). Moreover, the scaling functions governing 
the magnetic field and the temperature behavior of con- 
ductance fluctuations are found to be identical to their 



weak localization counterparts, see Sections III, IV 



The remainder of the paper is arranged as follows. Sec- 
tion [n] is devoted to the qualitative discussion of the role 
of the effect of the real electron-hole pair excitations on 
the weak localization and mesoscopic conductance fluc- 
tuations. The main point of this Section is to explain 
why the singlet excitations with transmitted frequency 
smaller than affect neither weak localization nor 

mesoscopic fluctuations. In Section HI we explicitly cal- 
culate the effect of interactions on mesoscopic fluctua- 
tions of conductance in one dimension, using the same 
approach as AAKLl We will also identify the diagram- 
matic contributions which are missed in the arguments 
for the role of r e <C in the conductance fluctuations. 
Section IV generalizes the calculation to two dimensions. 
Our findings are summarized in Conclusions. 



II. QUALITATIVE DISCUSSION 

The purpose of this Section is to explain interference 
processes, taking into account possibility of excitations 
of real electron-hole pairs, see also Ref. |l0|. For the weak 
localization correction, similar arguments were used in 
Ref. 0. 

A qualitative physical interpretation of quantum cor- 
rections is usually based on the following arguments, see 
e.g. Ref. O. Consider an electron diffusing in a good 
conductor, ppl 3> %■ Probability w for the electron to 
reach, say, point i starting from point /, see Fig. 0a, 0a, 
can be obtained by first finding the semiclassical ampli- 
tudes A a for different paths connecting the points, and 
then, calculating the absolute value of their sum, 



to 



(6) 



The first term in Eq. (0) is nothing but the sum of the 
classical probabilities of the different paths, and it may 



be found from the classical Boltzmann equation. The sec- 
ond term is the quantum mechanical interference of the 
different paths. In what follows, we will discuss the con- 
tribution of this term to transport and how it is affected 
by the electron-electron interaction. 



A. Weak localization correction 



For generic pairs a, j3, the product A a A% oscillates as 
the function of impurity configurations, see Fig. 0a. This 
is because the lengths of paths a and j3 are substantially 
different. As the result, contribution of such paths is not 
relevant for disorder averaged quantities but contributes 
to the mesoscopic fluctuations of the conductance. 

There are pairs of paths, however, which preserve the 
same phase, with the change of the disorder configura- 
tion. An example of such paths is shown in Fig. 0a. 
These paths almost coincide everywhere except the loop 
segment BEB (see Fig. 0a) which is traversed by trajec- 
tories I and 2 in the opposite directions. In the absence of 
the magnetic field and spin-orbit interactions, the phases 
of the trajectories 1 and 2 are equal. Therefore, the con- 
tribution of these paths to the probability w becomes 

\A 1 +A 2 \ 2 = \A!\ 2 + \A 2 \ 2 + 2ReA 1 A* 2 =4\A 1 \ 2 , (7) 

i.e. twice larger than the classical probability. Thus, in 
order to evaluate the weak localization correction to the 
conductivity, one has to determine the classical probabil- 
ity to find such a self-intersecting trajectory. 

Let us now consider the main effect of electron-electron 
interactions on the weak localization — excitation of soft 
electron-hole pairs. We consider processes involving ei- 
ther one excitation (probability Pi) or no excitations 
(probability Pq = 1 — Pi), see Fig. Allowing for 
the excitation of an electron-hole pair, one obtains 



A a — > A a + A . 



(8) 



where the superscripts and I correspond to the ampli- 
tudes involving emission of no electron-hole pairs or one 
electron-hole pair respectively. 

Because the states with different number of excita- 
tions are orthogonal to each other, we obtain, instead 
ofEq. @, 
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where the last two terms correspond to the interference 
correction. It is important to emphasize that the interfer- 
ence persists even if the final state contains an electron- 
hole excitation (last term). 
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FIG. 1. a) Example of classical paths between points i and 
/ contributing to the weak localization; b) The same paths 
with the excitation of one electron-hole pair with energy lo. 

We now notice that the emission of a soft electron-hole 
pair does not alter the geometrical form of the trajectory, 
thus, it does not change the classical probability corre- 
sponding to path a. As the result, we have 
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where amplitudes without superscript correspond to 
those in the absence of the interaction. What the emis- 
sion of the electron-hole pair may change, however, is the 
phase of the quantum amplitude. 

Indeed, denote the point of emission of electron-hole 
pair of energy lo on a classical trajectory by - time it 
takes for the electron moving along the trajectory a with 
energy ep to reach the emission point, see Fig. [I]. Denote 
the total time along the path a as t a . Than, the electron 
moves time t^ 71 with the energy e and time t a — t^ 71 with 
the energy e — lo. As the result, the geometrical phase, 
accumulated by electron, changes as 



argA^argAO-wCta-C 1 ) 



Thus, 
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Substituting Eqs. @ and into Eq. @, we obtain 
[instead of Eq. ([?])] for paths contributing to the weak lo- 
calization correction, 

\Al + A\+A% + A\\ 2 
= 2 | Ax | 2 + 2 | Ax | 2 [P + Pi cosw(i? m - t e 2 m )} ■ (12) 



The last term in Eq. ( |l2| ) describes the effect of the ex- 
citation of an electron- hole pair in the system on the weak 
localization correction. One can readily see that not each 
inelastic process destroys the interference. For instance, 
for lo — > 0, Eq. ( |T^ ) reproduces Eq. (^) exactly!! On the 
other hand, the time t^ m is shorter than r^. Thus, we 
may conclude that inelastic processes with energy trans- 
fer lo < I/t^ do not destroy the interference, which gives 
the physical reason for the low energy cut-off lo* ~ 1/t0 
in Eq. (@). 



B. Mesoscopic conductance fluctuations. 

Effect of inelastic processes. The arguments of the pre- 
vious subsection are easily generalized for the effect of in- 
elastic processes on mesoscopic conductance fluctuations. 
We can still talk about a pair of two paths, but now we 
will take those paths to be generic, see Fig. g. The inter- 
ference contribution from those paths, 



6G ~ 2ReAiAn 



(13) 



does not affect the average conductance because of ran- 
dom phases of those amplitudes, but it gives rise to the 
mesoscopic fluctuations of the conductance, 



(<5G 2 )^2(|A 1 | 2 )(|A 2 | 2 ). 



(14) 



Let us now consider the effect of the excitation of an 
electron- hole pair of energy lo. To do so, we use the qual- 
itative argument of previous subsection [sta rting from 
Eq. (|)] and substitute Eq. (|ll|) into Eq. (jl|). It yields 

SG ~ 2Re [AiA^(P + p ie ^(t2-ti-tr+*D)j . (15) 

Once again, we arrive to the conclusion that the ex- 
citations of frequencies smaller than the inverse times 
to traverse the trajectories, l/ii,2i do not change the 
interference correction. Similarly to the weak localiza- 
tion the lengths of paths are limited by r^. Thus, we 
may conclude that inelastic processes with energy trans- 
fer lo < I/t^ do not affect mesoscopic fluctuations, 
which gives the physical reason for the low energy cut- 
off lo* ~ I/t^ in Eq. (Q). Thus, inelastic time entering 
the weak localization and mesoscopic fluctuations should 
be approximately the same. The exact equality of those 
times will be proven in the next Section by a direct cal- 
culation, however, this result is definitely model depen- 
dent. Namely, it implies that the contribution of the 
quasi-static fluctuations in the systems does not over- 
whelm the role of the inelastic processes, and we discuss 
such fluctuations now. 

Effect of quasi-static fluctuations. In the linear re- 
sponse theory, a many-body system in its stationary state 
is excited at some time t\ and than the behavior of some 
observable quantity is studied at times t > t\. If the tem- 
perature is finite, the initial stationary state of the system 



3 



can be not only its ground state Eq, but also any of many- 
body eigenstates, E a \ the probability that the system is 
initially in such a state is oc e~ E "/ T . If there were no in- 
teraction, it would result only in the thermal average of 
the mesoscopic fluctuations. However, electron-electron 
interaction leads to the effective dependence of the disor- 
dered potential for electrons. The simplest, and the most 
effective example of this mechanism is the dependence 
of the Hartree potential of the electrons on the electron 
configuration. Since the measurable conductance is the 
result of the large number of measurements, each time 
the initial state may be different. 



a) 




b) 



sions, and u> <§; is the high-energy cut-off above 

which fluctuations can not be considered as quasistatic. 

To have a strong effect on mesoscopic fluctuations, the 
electric field should change significantly the wavefunc- 
tions of the electipns on the scale L<f,, which translates 
into the conditional 



e\SE\L^>^- = -. 



h 



(17) 



On the other hand, we estimate from Eq. ( |16| 



e*\6E(L<,)\'Li = TQ 



e 2 L 2 ~ d 



where g(L) is the dimensionless conductance on the lin- 
ear scale L. Taking into account Eq. (||) and the con- 
dition uiTfj, <C 1 we conclude that for the dephasing by 
the Coulomb interaction the condition ( JT7| ) can be never 
satisfied, and therefore the quasistatic fluctuations are 
negligible in comparison with the inelastic processes. We 
reiterate that this result does not hold for the scattering 
on the collective modes, which have peak in their spectral 
density on frequencies much smaller than I/t^. 
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FIG. 2. a) Example of classical paths between points i and 
/ contributing to the mesoscopic conductance fluctuations; b) 
The same paths with the excitation of one electron-hole pair 
with energy w. 

In principle, one could expect that the averaging over 
different configurations of the self-consistent potential 
may lead to an effect stronger than the excitations of 
the electron-hole pairs. This is possible, when there is 
an additional slow degree ofpfreedom such as magnetic 
impuritiesEa, moving defectsEJ, or slow fluctuations of 
the gauge fieldli-l However, this is not the case for the 
Coulomb electron-electron interaction, as we explain be- 
low. 

To find the magnitude of the effect, we first have to 
estimate the characteristic value of possible fluctuations, 
then evaluate the effect of such fluctuations on the meso- 
scopic fluctuations of conductance, and then compare 
this effect with effect of coming from the inelastic pro- 
cesses. According to Nyquist noise formula, the ampli- 
tude of the electric field 8E(L) fluctuating on the spatial 
scale L is given by 



5E 2 (L) 



T LU 



(16) 



where <7d is the conductivity of the system in d dimen 



III. CONDUCTANCE FLUCTUATIONS IN 
QUASI-ONE-DIMENSIONAL SYSTEMS 

In this Section, we consider a quasi-one-dimensional 
wire of length L and the number of transverse channels 
N± . The static conductance of the wire G is expressed 
through the non-local conductivity a(xi,x 2 ) as follows, 



72 J dxidx 2 (T xx (xi, x 2 ), 



(18) 



where x\ and x 2 label the coordinates along the wire. 
To simplify the expressions, we disregard first inelastic 
processes and include them later on. We express the 
symmetric part of the conductivity in terms of— Glen's 
functions and substitute it in Eq. (18). We find' 16 



G 



dr\dr 2 



-f~3 Xl G R ( ril r 2 ;e)j X2 G A < 
7r de 



r2,rv,e 



(19) 



where the integration is performed over all the sample, 
the spin degeneracy is taken into account, / is the Fermi 
distribution function, and the current operator j x is de- 
fined as follows, 



913x92 



2m 



(92d x gi - g\d x g 2 ) 



For the rest of the article, we employ the system of units 
with h = 1, and restore h in the final results. 
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(a) 



(b) 



FIG. 3. Diagrams with two diffusons (a) and two cooper- 
ons (b) contributing to the conductance fluctuations. Dashed 
lines represent impurity scattering, interaction is not yet 
taken into account. 

In the following, we consider only high-temperature 
limit, T 3> D/L 2 , because it is the only case when the 
inelastic processes (rather than sample size, L) are con- 
trolling the magnitude of the fluctuations. In this case 
the main contribution to the conductance fluctuations is 
given by two-diffuson and two-cooperon diagrams, Fig. [|. 
The resulting correlation function for the conductance 
fluctuations at different magnetic fields Hi, H2 is ex- 
pressed in the time domain as 



6G(Hi)5G(H 2 ) = 



(2e 2 D) 2 
SttTL 4 



dxidx2 I dt 



(20) 



\V%( Xl ,x 2 ,t)\ 2 + \V 1 c ?{xi,x 2 ;t)\ 2 



where the overbar stands for the disorder averaging. De- 
riving Eq. (Boh, one makes use of the approximation 



dei de 2 
2tT 2tt" 
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justified at time scale larger than 1/T. 

Semiclassical retarded diffuson and cooperon propaga- 
tors entering into Eq. p0|) are solutions of the equations 



d t - Dd 2 x 




X j X 1 

V£(x,x';t) 



S(x-x')S(t), 



(21) 



where D is the diffusion coefficient, and the symmetry 
breaking parameters T(? D are defined as (see Ref. |?]) 



1 _ e 2 a 2 D(Hi - H 2 ) 2 1 



e 2 a 2 D{Hi+H 2 f 



I2fi 2 c 2 



(22) 



with a being the transverse dimension of the sample. It 
is worth mentioning that the numerical coefficient here 
is geometry dependent. 

So far, we merely followed a standard avenue (see e.g. 
Ref. [r]). Now we are prepared to introduce electron- 
electron interactions. On the language of diagrams, we 



must add to Fig. || all of the possible interaction lines. 
Since inner and outer rings represent the measurement at 
significantly different times, the interaction lines do not 
connect these two rings, and only may be drawn within 
the same ring, connecting G R with G R , G R with G A , and 
G A with G A for the same impurity configuration. Fol- 
lowing Ref. [l], these lines are conveniently represented by 
external time-dependent random fields, ip a (x,t), where 
the index a assumes values a = 1 (outer ring) and a = 2 
(inner ring). These fields are assumed to be Gaussian 
distributed with zero average. The correlation function 
is described by the Keldysh component of the propagator 
of the screened Coulomb interaction, 



(<p a (x,t)^(x',t')) 



5 a p5{t - t') 



2T 



Dvi J 2?r q 2 



dq 1 iq(x-x') 



(23) 



where V\ is the thermodynamic density of states per unit 
length. Equation (|23| ) is nothing but a space-time version 
of 



Tm 



2T Dq 2 



Dq 2 vi 



and we assumed T > u. This assumption is justified, 
because the main contribution to the dephasing rate is 
coming from the energy transfer oj much smaller than T. 
(The diagrams explicitly showing cancellation of all the 
processes with u > T can be found, e.g., in Refs. p^|,|6[) 
Because we also disregard all effects due to finite size of 
the sample, this implies the following hierarchy of energy 
scales, 



T » tV 1 » E c = D/L 2 



(24) 



(here T<f> stands not only for the phase-relaxation time, 
but for all time scales due to electron-electron scatter- 
ing). In the following, we assume that the conditions 
( p4| ) are satisfied. 

The factor <5 Q/ g in the right-hand side of Eq. ((2^) explic- 
itly indicates that the fields attached to outer and inner 
rings of the diagram Fig. || are uncorrelated, i.e. no in- 
teraction lines, indeed, can be drawn between the rings. 
The momentum integral in Eq. ( p3| ) diverges, but our fi- 
nal result will contain well-defined differences of integrals 
of this type. 

Introduction of the fluctuating fields modifies the equa- 
tions for the diffuson and cooperon (|Tj), see Fig. ^ and 
Ref. [I], which now become the functionals of the fluctu- 
ating fields, 



d t -Ddl+i(ip a {x,t)-<pP{x,t)) + 



V^(x,x';t;{ip a (x,t),^(x,t)}) 
Vf(x,x';t;{^(x,t),^(x,t)}) 
= S(x-x')S(t). 



(25) 
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The correlation function of conductances is given by 
the equation similar to Eq. (|20|), but all the interaction 
lines in Eq. are connected by the propagator (23) 



SG(Hi)5G(H 2 ) = ( ^ T ^ 4 J dx x dx 2 J dt (26) 

(\V 1 D 2 (xux 2 ,t)\ 2 ) v + (\Vl?(x 1 ,x 2 ;t)\ 2 ), 

where (. . .) v stand for the averaging over the fluctuating 
field ip 1 ' 2 . 



R, a 



x,t 



x',0 



A, p 



CF diffuson V^f ■ Zigzag lines represent random 



FIG. 4. 
fields <p a >t 



Before we perform actual calculation in Eq. (|26j), we 
pause for a moment to discuss a-j"elation of this formula 
with the other theoretical workQ. We observe that the 
propagator (Pjf}^ contains all possible interaction lines 
drawn between G R and G R , and also between G A and 
G A , but not between G R and G A . This is exactly an 
object (let us call it CF diffuson), which determines the 
out-scattering term in the collision integral in the Boltz- 
mann equation, and it was studied in details in Ref. |6|. In 
contrast to the "ordinary" diffuson, which is insensitive 
to electron-electron interaction due to Ward's identity 
(charge conservation), the CF diffusion (T'}^) ip acquires 
a massive pole, real part of which is identified with the 
out-scattering time r e . One can thus imagine (and this 
was, indeed, conjectured in Ref. ^) that the temperature 
dependence of conductance fluctuations is governed by 
the time r e , which is parametrically different from r^. 
The calculation presented below shows that this conjec- 
ture is not correct. The resolution of this fallacy is that 
the averaging in Eq. (prj|), which is essentially coupling 
of all random fields (p a according to the rules (p3|), pro- 
duces not only a contribution which contains averages 
\{V 12 ) V \ 2 (Fig. §t), but also dia grams where interaction 
lines connect upper and lower Green's functions within 
the same ring (Fig. |B|b). Both contributions diverge in 
the infrared limit (and have to be regularized in order to 



extract sensible resultsQ), but their sum is well-behaved. 

To proceed with the evaluation of Eq. (p6|), we write 
V a P as a functional integralliao, 



V^ c (x,x';t;{(p(x,t)}) 



0(t) 



e °.o 



y(t)=x 



Vy(r) 



< \p [ J dr 



.&+i^[y(r),r}-i^[y(r),T^y 

(27) 



where 9(t) is the step function, Z is the normalization 
factor, that will be included in the measure of the func- 
tional integration in all of the subsequent formulas. Sub- 
stituting this expression into Eq. (^6|), and averaging over 
Gaussian random fields ((e lif ) — e~*> v ^ 2 ), we obtain 
with the help of (||) 

5G(ffi)<SG(ff 2 ) = V 3?rTL ; 4 J d Xl dx 2 dt (28) 

yi(t)=x 1 rV2(t)=x 1 

V Vl (t) / Vy 2 {t) 

yi(0)=X2 J V2(0)=X2 



x e o + e a 



x exp 



dt' 





(a) 



(b) 



FIG. 5. Examples of diagrams with interaction (shown as 
zigzag lines) contributing to conductance fluctuations. The 
diagram (a) is reduced to the CF diffuson, while the diagram 
(b) is not. Conclusion about the differences of inelastic rates 
for weak localization and conductance fluctuation is a conse- 
quence of missing the diagram (b). 

Following Ref. [I], we introduce new variables, 
yi(t)±y 2 (t) 



«i,a(*) = 



V2 



This yields 



(2e 2 D) 2 

8G{Hi)8G(H 2 ) = \ ^J A I dx x dx 2 



xJ 1 (V2x 1 ,V2x2;t)J 2 (0,0;t) [e^ + e" 7 ? 



(29) 



G 



nZi{t)=Xi 

J±(xi,x 2 ;t) = / Dz%(t)exp 

«/*l(0)=s 2 
rz- 1 {t)=x 1 

J 2 (x 1 ,x 2 ;t) = / Dz 2 (t) 

Jz 2 (0)=X 2 



1 7 2 



x exp ■ 



d£ 



4D 



2y/2Te 



where o\ is the one-dimensional conductivity, and we 
used Einstein relation a\ = e 2 viD. 

Now we represent these functional integrals Ji,2 as so- 
lutions of differential equations. The integral J\ solves 



(d t -Dd^) J 1 = S{t)S(x 1 -x 2 ). 



(30) 



Integrating both sides of Eq. ( |30| ) over x\ and x 2 , and 
neglecting the boundary term at t <C L 2 /D, we obtain 

J dx 1 dx 2 Ji(V2x 1 ,V2x 2 ,t) = -j=9(t). (31) 

Similarly, J 2 obeys the equation 



fit - Ddl + 



2 v / 2Te 



-\x 1 \\j 2 = 5{t)5{x l -X2). (32) 



Substituting Eq. @ into Eq. @, and using Eq. (|3|) 
we find 

8G(Hi)5G(H 2 ) = = 0) + Q c (:r = 0)) . 



and Qd,c{x) obeys the equation 



(33) 



r 12 
'C.Z3 



-Ddl 



2 v / 2Te 

0"! 



■N Qo, D (aj)=*(i). (34) 



Equation has been previously considered in Ref. [j], 
and it has the solution in terms of the Airy function 

Ai(x), 



,d(x) 



Ld 



A 



\ T C D 



\/2\x\ 

DLa, 



2V2D AV (_^_) 



where the dephasing time and the dephasing length 
have exactly ike same form as for the weak local- 
ization correctioni'B (numerical coefficient is corrected in 
Ref. 0), 




(36) 



Substituting Eq. (|35|) into Eq. (|33j), one finally obtains 



5G(H 1 )5G{H 2 ) = { - 



hD L 
3irL 2 T~L 



T 12 

T D 



It. 

T 12 

~c 



T)(x) = - 



1 



[]nAi(x)}' 



(37) 



Equation ([57|) with entries (|2^) and ( pq ) is the main 
quantitative result of the present Section. It shows that 
the dephasing rate governing temperature and magnetic 
field dependence of the mesoscopic fluctuations is exactly 
the same as in weak localization. Moreover, this result 
can be combined with the expression for the weak local- 
ization correction 



SGwl(Hi) 



5a W L (Hi) 



6a WL (Hi) 



irk 



->1 



It. 

Tq 



to the form free of geometrical uncertainties (as well as 
uncertainties in the value of the diffusion coefficient) , 



sa ( iio<mii^ = r T )^ 



(38) 



SG 



WL 



Hi — H 2 



SG 



WL 



Hi + H 2 



This result gives the relation between two measurable 
quantities, and thus may serve as a test for the dephasing 
mechanism. Equations (p^) and (38) are valid provided 



^/ T L)C K is also assumes that there is no spin-orbit 

interaction. It may be shown that in the case of strong 
spin-orbit (SO) interaction, the result (|3|) still holds up 
to a numerical factor of 1/2. In the case of the crossover 
between strong and weak SO interaction one has to iden- 
tify the singlet 5G S and triplet SGt contributions to the 
weak localization correction SGwl = 3SGt — SG S by 
corresponding fits and replace SGwl in Eq. ( |38| ) with 
[SG WL + 2SG S ] /2 = [3SG t + SG S ] /2. 

Now, for conceptual clarity, we employ the result (J37 
to extract the relaxation time associated with conduc- 
tance fluctuations. It is important that this time is un- 
physical by itself, and only has a meaning when explicitly 
linked to Eq. @. 

For this purpose, we take Hi — H 2 — and define the 
time tt as a mass in the pole in the CF diffuson and 



(35) CF cooperon V^ 2 which enter Eq. (|fj). Writing 



V 1 c 2 n(x,x',t) = 



(27T)* 



Dq 2 



iu> - 



substituting this expression into Eq. (|20j) and performing 
the integration, we obtain for conductance fluctuations 



sc* ^ 2 

6TTHTL 3 



(HL 
\D 



(39) 



Comparing this to the result (37), we identify the inelas- 
tic relaxation time tt responsible for the temperature 
dependence of conductance fluctuations, 
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t t =?7 2 (O)t «O.53t 



(40) where the vector potentials are such that 



where is defined in Eq. (|3q), i.e. it is precisely the 
same time one obtains if one considers weak localization 
by introducing a finite mass in the pole of the Cooperon. 
Thus, the temperature dependence of conductance fluc- 
tuations does not produce a new time scale as compared 
to Eq. ( |36| ) and is certainly not determined by the out- 
scattering time r e . The numerical coefficient 0.53 reflects 
the behavior of the scaling function ( p7| ) in low magnetic 
fields. 



IV. TWO-DIMENSIONAL CASE 

Equation (|3^) can be readily generalized to the two di- 
mensional sample, and we outline the main steps of the 
corresponding derivation. 

Consider a two-dimensional system of the size L. Per- 
forming the same steps as in the derivation of Eq. ( |26| ) 
one finds: 



8G{H X )8G[H 2 ) 



{2e 2 Df 



x [(\V^{n,r 2 ,t) 



3ttTL 4 

2 S 



If 



d 2 r x d 2 r 2 
rh 2 {r 1: r 2 ;t) 



dt 

1 2. 



(41) 



where two-dimensional integrations are performed within 
the sample, (. . .) v stand for the averaging over the fluc- 
tuating field ip 1,2 with correlation function analogous to 
Eq. A 



(^(r,iW,t')) 

2T f d 2 q e iq( - r 

8 a p6(t-t')- 



Du 2 J (2tt) 2 



(42) 



with v 2 being the thermodynamic density of states per 
unit area. In Eq. (|42[), the integration is limited from 
above by |g| ~ (T/Dy/ 2 . Such an accuracy of the ultra- 
violet cut-off is sufficient for the logarithmically divergent 
integral. 

Diffuson and cooperon propagators entering Eq. jil] ) 
are the solutions of the two-dimensional analog of 
Eq. (§), 



dt-D 



7 q/3 



+i(p«(r,t)-^(r,t)) + 



V^(r,r';t;{^(r,t),^(r,t)}) 
= 8(r-r')6(t), (43) 

where times 1/td,c, see Eq. (p2|), describe the effect of 
the magnetic field component parallel to the film plane. 
The effect of the magnetic field perpendicular to the 
plane is described by 



v a/3 = 
v 7 — 

ACtfi _ 

D — 



a, = 1,2; j = D,C; 



A a - A 13 ; 



A^f = A° 



A? 



(44) 



V x A a = , 



and Hj_ is the component of is the magnetic field per- 
pendicular to the plane. 

Transformations leading to Eqs. (^6|) and ( p5|) are 
pretty much the same as in ID provided we make ob- 
vious changes x — > r, q — > q, d x — > V r . Writing again 
the CF diffusons and cooperons V%f c as functional in- 
tegrals (|27| ) and performing an averaging over Gaussian 
fields (p a , we obtain a two-dimensional analog of Eq. (pE}), 

(Op 2 ]!) 2 r r°° 

8G(Hi)8G(H 2 ) = ^-^jt J d 2 r x d 2 r 2 J dt (45) 



_ 2t rVi(t)=ri rV2(t)=ri 

x J2 e w \ v v^) l v v*Q) 

~t=D,C •'Vi(0)=ra 
rt 

dt' 

4D 4D 



V2(0)=r 2 



x exp 
ie 



• 9 "2 



A 12 { yi ) yi -A 12 {y 2 )y 2 ) 



4T 



dq 1 



Dv J (2tt) 2 q 2 
Introducing new variables 

R( t) = m+Mi 



[1 - co S (q( yi (t') - y 2 (t'))} 



r(t)=yi{t)-y 2 (t), 



and reducing the functional integrals back to differential 
equations, we obtain the result 



6G{H{]SG{H 2 ) 
_ (2e 2 D) 2 



[Q D (\r\=L T ) + Q c (\r\=L T )}, (46) 



SttTL 2 

where Qd,c{ t ) obeys the equation 

1 



-D(V 12 f + U(r) 



T 12 

T D,C 



and the potential is given by 

2T f dq 1 - cos(qr) 



Dv J (2tt) 2 



Qc, D (r) = 8(r), (47) 



T , [L T + r 
■ In 



ttDv 



(48) 



where the last expression and Eq. ( }46| ) are written with 
the logarithmic accuracy and we take into account the 
high-momentum cut-off at q ~ L^, , Lt = (D/T) 1 / 2 . 

Equations ( fltf ) - ( [iq ) should be compared with the 
corresponding expressionpfftr the weak localization cor- 
rection in two dimcnsionso'El, 
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Sa(H 1 



C(r = l), 



(49) 



-D(V 1 c 1 ) 2 + U(r) + ^ T 
T c 



C(r) = 5(r), 



where the logarithmic divergence should be cut at the 
elastic mean free path I. 

Therefore, we conclude that the relation similar to 
Eq. (|38|) should hold, 



SG{H = 0)5G{H = Q) - 5G(Hi)8G(H 2 ) 



(50) 



8 a 



hD 

h J 3L 2 T X 
f Hi — B.2 



WL 



5°wl ( Hl ~l Hl ) - 2SawL (0) 



It is important to emphasize that the relation (|5C 
holds even before one starts an approximate solution of 
Eq. (^). Note however th at the r esult similar to Eq. (H|) 
does not hold, since both 5GSG and SawL diverge log- 
arithmically with different cut-offs. This is why in Eq. 
( pp| ) we had to subtract zero-field contributions, which 
cancels logarithmic divergences. 

The effect of the spin orbit interactions on our final 
result ( |50| ) is the same as for one-dimensional geometry, 
see discussion after Eq. (38). |-, n 

We write here the explicit expressionlirQ for the weak 
localization correction in two dimensions for the reference 
purpose, 



6a WL (H,T) 



2ir 2 h 



In ■ 



1 



H 





- 1 )} 







where ^(x) is the digamma function, £Ih — 4:eDH±/ch, 
and r* is determined by the equation 



f 



1 

th 



Te 2 Rr 



■ In 



T 



h 2-kTi h/r* + hn H ' 



Similarly to one dimension, we can also extract the in- 
elastic time Tt, defined as a pole of CF diffuson in zero 
magnetic field. An explicit calculation gives tt ~ t^. 
This relation contains a numerical coefficient of order 
one, which can only be determined by going beyond the 
logarithmic accuracy. We do not attempt such a calcu- 
lation in this paper. 



comparing it with the weak localization data obtained on 
the same sample may give information on the nature of 
inelastic interactions in disordered metals. 

We are not aware of attempts to make such a compar- 
ison between inelastic times directly. However, recently 
Hoadley, McConville, and Birge (HMB)E3E3 presented 
very careful measurements of the magnetic field depen- 
dence of l//-noise in silver films. A standard assump- 
tion in the theory of 1/ /-noise in metals (for review, see 
Ref. |2l|) is that it is produced by low-frequency motion of 
impurities. Mathematically, the magnitude of 1/ /-noise 
in such a model is given by a set of diagrams identical 
to those for conductance fluctuations (Figs. ||, ||) with 
the only difference that external and internal -rings are 
described by different impurity configurationsE3c3. As 
the result the field dependence and the temperature de- 
pendence of the noise should be given by the paramet- 
ric derivative of Eq. (|50|), i.e. it should be expressed 
through the derivatives of the parallel field dependence 
of the weak localization. 

HMB compared the timescale defined as a pole in the 
diffuson (in our notations, Tt), with the phase relaxation 
time 70, extracted from their own measurements of the 
weak localization correction on the same films. Their 
procedure results in tt — which was interpreted 

to be consistent with the theory of Ref. 0. Our results 
(|50| ) contradict that interpretation. 

To our opinion, the only possible reason of this dis- 
crepancy is the electron-electron interaction in the triplet 
channel which we did not take into account. This inter- 
action can be singled out in experiments with the materi- 
als with stronger spin-orbit scattering. Other sources of 
l//-noise seem to be excluded, since the functional form 
of the experimentally measured by HMB magnetic field 
dependence perfectly fits theoretical predictions. Dc- 
phasing on slow moving impurities itself, see discussion 
in Sec. II B, would give a temperature dependence differ- 
ent than that in experiment and may be ruled out. Wc 
believe that the contradiction between the theory and 
the experiment revealed in our paper indicates that the 
quantitative study of inelastic processes in mesoscopic 
samples remains an interesting topic and deserves future 
investigation. 
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